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ABSTRACT 

The study of extragalactic integrated light can yield partial information on stellar 
population ages, abundances, and the initial mass function (IMF). The power-law slope 
of the IMF has been studied in recent investigations with gravity-sensitive spectral 
indicators that hopefully measure the ratio between KM dwarfs and giants. We explore 
two additional effects that might mimic the effects of the IMF slope in integrated light, 
the low mass cutoff (LMCO) and a variable contribution of light from the asymptotic 
giant branch (AGB). We show that the spectral effects of these three (IMF slope, 
LMCO, AGB strength) are very similar and very subtle compared to age-abundance 
effects. We illustrate parameter degeneracies and covariances and conclude that the 
three effects can be disentangled, but only in the regime of very accurate observations, 
and most effectively when photometry is combined with spectroscopy. 

Key words: galaxies: abundances — galaxies: evolution — galaxies: elliptical and 
lenticular, cD — galaxies: luminosity function, mass function 


1 INTRODUCTION 


Deriving astrophysical parameters from integrated-light ob¬ 
servables is a widely accepted practice in the field of 
Galactic and extragalactic research. However, increasing 
evidence shows that single-burst, single-composition stel¬ 
lar p opulations oversimp l ify the underlying stellar sys¬ 
tems JCratton et al.l l20ld : iKavirai et al.l l2007j l . Additional 


parameters in the modelling process, such as: multiple 
burst-age stellar popu lations (ITrager et al.l l2000al . l2005l : 
iGoudfrooii et all 12011). metahicity and helium abundance 
variation lNorri'gl2004l : |Lee_et_|^ 2005ll. init i al mass func¬ 
tion (IMF) va riation ( Weidner et al. I I2OI3I: iBekkil I2OI3I: 
|ChabrieL_eL_^ 20.1 4| ), and chemica l abundance variation 
I Potter et al.l l2007l : (bee et al.l l2009l l were investigated in 
hopes of reconciling various contradictions between ob¬ 
servation and theory. This group explored metallicity- 
compositeness by varying the abundance distribution func¬ 
tion (ADF; dM/d[M/H]; the mass fraction of the stellar 
population at each [M/H]) effect in iTang et al.l (2014, Pa¬ 
per I). The ADF shape is a gentle rise at low abundance, a 
peak, and a steeper fall-off at high abundance. By varying 
the widths of the ADFs we discovered “red lean” and “red 
spread” phenomena. “Red lean” means that a narrower ADF 
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appears more metal-rich than a wide one, and “red spread” 
describes that the spectral difference between wide and nar¬ 
row ADFs increases as the ADF peak is moved to more 
metal-rich values. 


We continue to explore in this paper additional un¬ 
derexplored effects, this time ones that might endanger 
measurements of the IMF slope in integrated stellar pop¬ 
ulation (SP) models. The IMF is of importance in stel¬ 
lar and extragalactic astronomy, affecting galaxy luminos¬ 
ity evolution, star formation, chemical evolution, and mass 
budget between normal and dark matter. The IMF regu¬ 
lates the mass distribution of the stellar population, and 
thus impacts the luminosity function, mass to light (M/L) 
ratio, and number of stellar remnants. However, direct 
IMF slope derivations that counts individual stars are lim¬ 
ited by observational capabilities, uncertainties concerning 
the mass-luminosity relation, stellar evolution, dynamical 
evolution, binary fraction, and many other factors. Re¬ 
cently reemerged gravity-sensitive spectral lines (e.g., GaT, 
Nal, and Wing-Ford band) show bright future applica¬ 
tions for indirectly estimating the IMF slopes of unresolved 


galaxies dGenarro et al. 

I2OO3I: van Dokkum & Gonrovll2010l: 

iGonrov & van Dokkum 

l2012bl). 

Gravity-sensitive spectral 


lines measure the dwarf/giant ratio. Our stance is that this 
ratio could be altered by effects other than the IMF slope. 
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2 Baitian Tang and Guy Worthey 


First, a low-mass cut-off for the IMF (LMCOj^i exists 
in every set of SP models. It is the mass limit of stars in¬ 
cluded at the lower mass end — any star with mass smaller 
than this limit is assumed to have no contribution in the SP 
models. Often, the LMCO is a pragmatic choice dictated 
by the choice of stellar evolutionary isochrones that go into 
the SP models. It is readily appreciated that the LMCO is 
closely related to the fraction of low mass stars, and thus 
the dwarf/giant ratio. Therefore, we may encounter degen¬ 
eracy when determining IMF slope and LMCO simultane¬ 
ously. The LMCO is variable among different SP models 
found in the literature. Fo r example, the Padova isochrones 
dBertelli et ahllioosl . l200^ set 0.15 Mp, as the LMCO, while 
the composite isochrones of IConrov fc van Dokkuriil (l2012al l 
sets 0.08 Mq as the LMCO. As a result, the derived IMF 
slope values may change when swapping between SP models. 

Second, we note the possibility that the dwarf/giant 
ratio may come as easily from modulating the number of 
giants as modulating the number of dwarfs. The bolometri- 
cally brightest giants are asymptotic giant branch (ACB) 
stars. These evolved stars have long be en known to be 
diffic u lt to constrain in most SP models llConrov fc Gumil 
I 2 OI 0 I : ICirardi et al.l I 2 OI 0 I . I 2 OI 3 I I due to uncertain stellar 
evolution (in turn due to uncertain rules about mass loss 
in giants) and also di fficult to empirica ll y constrain due 
to counting statistics dFrogel e t al.l1l99Cll: ISant os fc Fro^ 


ll997l : [Bruzual fc Charlotll2003l : ISalaris et al.ll20l4l . We sug¬ 

gest a plausible, hypothetical trend, ADF-AGB-IMF mas¬ 
querading, in which increased metallicity causes cooler gi¬ 
ants, which causes increased mass-loss, which causes fewer 
AGB stars, which resembles an increase in IMF slope. 

In this paper, we study spectral and photometric vari¬ 
ations in response to simultaneous changes of IMF slope, 
LMCO, and AGB strength (321 • We show that the degenera¬ 
cies can be marginally lifted for old, metal-rich populations, 
but the degeneracies are still firm in young populations ( H2.2\ 
& tl2.3l) . Next, we dig into the mechanisms behind all these 
variations by studying the numbers of stars, luminosities, 
colors, and spectral indices of each evolutionary phase ( W2A\i . 

Bottom-heavy IMFs are indicated i n local mas¬ 
sive and metal-rich elliptical galaxies (iGenarro et al 


200 ^ yaB_Dokkum_&_Cpnro 2 2010 l : [^nrov fc van Dokkum 

2012b r TCappellari et al.l 12012 ). On th e other hand , star- 


formation theorie s such as iLarso 3 ll 19981. I 2005 II and 
iMarks et al.l (l2012l ~l im ply the metal-fr e e earl y Universe fa¬ 
vors a top-heavy IMF. IWeidner et al.l (l2013l l points out a 
time-independent bottom-heavy IMF generates too little 
metal and fewer stellar remnants than observed. The hy¬ 
pothesis that IMF steepens as the Universe evolves will 
change the light ratio between metal-poor and metal-rich 
stars in a similar way as the ADF effect, thus we call here 
ADF-IMF coupling. In 33 we build rough models with ADF- 
IMF coupling to explore the ramifications, at least qualita¬ 
tively. The coupled models appear more metal-rich than the 
noncoupled models, due to the suppression of metal-poor 
stars, resembling the red lean effect. 

Finally, we test parameter recovery using a Monte Carlo 
approach. The mean recovered values agree with the input 
values inside the error range without alarming systemat- 


ics ( 1)4.41) . We uncover covariances among the IMF slope, 
LMCO, and AGB effects when combined with the more 
usual age and metallicity effects. Though the magnitudes 
of the IMF-related effects are smaller than the latter ef¬ 
fects, these two groups of effects vector almost orthogonally 
(323J, and we prognosticate bright hopes for disentangling 
all of these stellar population parameters. A summary of our 
results are given in 33 


2 IMF SLOPE, LMCO, AND AGB 
2.1 Model description 

A new version of o ld integrated-light models (IWorthevlll994l : 
iTraeer et al.l[l998h is adopted. The new models use a n ew 
grid of synthetic spectra in the optical (iLee et al.ll2009l l in 
order to address the effects of changing the detailed ele¬ 
mental composition on an integrated spectrum. The models 
retain single burst age and metallicity as parameters, but 
were expanded to also include metallicity-composite popu¬ 
lations, and the three IMF-related parameters we discuss 
(IMF slope, LMGO, and AGB modulation). 

F o r this work, we adopt the isochrones of iBertelli et al.l 
ll2008l. l2009l'l using the thermally-pulsing a symptotic gi¬ 
ant br anch (TP-AGB) tr eatment desc r ibed i n iMarigo et al.l 
ll2008ll . Improving upon IPoole et al.l ll2010ll . stellar index 
fitting functions were generated fro m indices m e asured 
from the stellar spectral libraries of IValdes et al.l (l2004fl 
as re- fluxed by 1. Chilingar ian (2014, private communica¬ 
tion) , IWorthev et al. ll2 014ah. M ILES dFalcon-Barroso et al.l 
l201lil . and IRTF ( Ravner et al.l[2009ll . all transformed to a 
common 200 km s“^ spectral resolution via smoothing if the 
spectral resolution was greater or by linear transformation 
after index measurement if the resolution was lesser. Index 
errors were assigned by iteratively comparing the libraries 
and seeking error models of the form = <raoov + C'^photom 


where Ufloor is an error floor mostly applicable to indices of 
large wavelength span, and C is a fitting constant. Multi¬ 
variate polynomial fitting was done in five overlapping tem¬ 
perature swaths as a function of O^ff ~ 5040/Te//, log g, 
and [Fe/H]. The fits were combi ned into a lookup table for 
final use. As in IWorthe'J (ll994lL an index was looked up 
for each bin in the isochrone and decomposed into “index” 
and “continuum” fluxes, which added, then re-formed into 
an index representing the final, integrated value after the 
summation. 

The IMF slope, LMCO, and AGB effects are calculated 
at two representative values. The IMF slopf|3 is calculated 
at 2.35 and 1.70, the LMCO at 0.15 and 0.30 Mq, and the 
AGB contribution with full and 80% strength. For model¬ 
ing, we define the AGB phase as beginning 0.5 mag brighter 
than the red clump, and continuing to the end of the star’s 
life. The number of AGB stars of each mass bin is reduced 
to 80% to simulate a weaker AGB strength model. The SP 
model with IMF slope of 2.35, LMCO of 0.15 Mq, and full 
AGB strength is chosen as the standard model. Parameter 
variations from these nominal values are carried with lin¬ 
ear interpolation. Naturally, if the parameters drift too far 


^ Throughout this paper, IMF LMCO is called LMCO for short. 


^ A power-law function, where 2.35 is the ISalpeterl lll955ll slope 
and 1.70 is a bottom-light hypothesis. 
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On Disentangling IMF Degeneracies in Integrated Light 3 


from the modeled values they can no longer be considered 
realistic. 

To illustrate the different effects we vary one parameter 
at a time. That is, the results we display graphically are 
partial derivatives (dindexldeffect)- 

2.2 An Old, Metal-rich Population 

Since most of the elliptical galaxies in the local Universe are 
old and metal-rich, we pick age = 10 Gy, and \og{Z) — 0.37 
to mimic a typical elliptical galaxy. In the left column of 
Figure [T] we show optical-near infrared (NIR) color-color 
plots. Note that we set the same y axis dynamic range (0.1 
mag) for all the color-color plotqj. An extinction vector of 
Av = 0.02 mag is sketched on the bottom right of each 
plots. We note that (1) the IMF slope, LMCO, and AGE 
effects are not parallel vectors, which means it is practical to 
lift the degeneracy of a single index with a mix of different 
colors. However, we note that these effects are detectable 
at the level of ~0.02 mag. Such accuracies are comparable 
to contemporary detection limits and therefore technically 
challenging by today’s standards. In addition, (2) the {J—K) 
vs. (B — V) plot shows the weakest drifts among the three 
plots, due to the wavelength-proximity of these two pairs of 
filters. This wavelength-proximity is also the reason for small 
extinction vector. We also find that (3) decreasing the IMF 
slope and AGE strength tend to change the {X — K) colors in 
the opposite direction. Interestingly, increasing the LMCO 
leads to a decrease in {V — K) and {R — K), but an increase 
in (J — K), implying blue (J — K) colors for stars between 
0.15 and 0.30 Mq f i)2.4ll . (4) Finally, using colors alone, dust 
extinction is almost indistinguishable from LMCO effects. 

In the middle and right columns of Figure [T] we plot 
[MgFeQ versus five IMF-sensitive indices. We introduce to 
the Worthey m odels the optica l IMF -sensitive indices bTiO 
and CaHl from lSpiniello et al.l (|2014| L and plan to use them 
for future intermediate-redshift IMF research (Tang et al., 
in preparation). We notice that (1) the drifts caused by 
weakening the AGE strength are the smallest, due to low 
flux contribution of the AGE stars at the age of 10 Gyr 
mM- We also see that (2) decreasing the IMF slope tends 
to decrease the [MgFe] index (about 0.02 A), while increas¬ 
ing the LMCO leads to small variation of [MgFe]. In ad¬ 
dition, (3) the smaller dwarf/giant ratios induced by the 
IMF slope and LMCO effects lead to smaller values for IMF- 
sensitive indices. Eut the magnitude ratios of these two ef¬ 
fects vary for different indices. The results that the LMCO 
drifts are larger than t he IMF slope dr ifts in Ti02 and WFE 
A (Wing-Ford band) (Wing fc FordI 1 19691 : IWhitforJ 1 19771 : 
iHardv fc Coutura Il98f ) suggest caution should be taken 
when interpreting these two indices as a clean IMF slope 
indicator. 


2.3 A Young, Metal-rich Population 

ISchiavon et al.l ll2006ll suggested young elliptical galaxies 


® 0.01 mag and 0.15 A for magnitude and angstrom unit indices, 
respectively. 

4 [MgFe] = ^Mgb * (Fe5370 + Fe5335)/2) 


with ages of the order of 1 Gyr prevail in the early Uni¬ 
verse (z~0.9). Studying these galaxies places important con¬ 
straints on the elliptical evolution models. Here we pick age 
= 2 Gyr, and log{Z) = 0.37 to mimic a young elliptical 
galaxy. In the all the plots of Figure [2] weakening the AGE 
strength induces a stronger drift than the old population, 
due to the hig h flux contribution of AGE stars at the age 
of 0.2 — 2 Gyr llMaraston et al.ll2006l L To accommodate the 
large AGE drifts, the y axis dynamic ranges are set to 0.25 
mag for color-color plots. Note that this may causes visual 
differences between the color-colors plots of Figures [T] and (2] 

In the optical-NIR color-color plots of Figure [21 (1) in¬ 
creasing the LMCO leads to negligible signal, while weaken¬ 
ing the AGE strength shows a drift that is a factor of two 
larger as the old population. Another noticeable change is 
(2) the optical colors on x axis. Instead of getting bluer as 
the old population does, decreasing the IMF slope of a young 
population drives the [B — V) and [B — R) colors redder, 
probably due to stronger post main sequence phases f 1)2.41) . 
Comparing with the extinction vectors, (3) the AGE drift 
and IMF slope drift should be clearly detectable at 0.02 mag 
level. 

In the middle and right columns, (1) the index increases 
of bTiO, Ti02, and WFE A caused by decreasing the IMF 
slope of a young population confirm that IMF determination 
is sensitive to the age parameter. We also see that (2) the 
shallower IMF slope model also has greater [MgFe] than the 
standard model. In addition, (3) the AGE and IMF slope ef¬ 
fects vector almost oppositely in most of the Figure (21 plots. 
The partial derivative nature of our results reveals the de¬ 
generacy of IMF slope and AGE strength is still firm for the 
young population. We further discuss the reasons for dif¬ 
ferent drift directions by inspecting the luminosity weights, 
colors, and indices of each evolutionary phase in 1)2.41 

2.4 Colors and Indices Broken into Evolutionary 
Phases 

To illustrate the reasons behind the color drifts of three ef¬ 
fects, we show the number of stars, luminosities, and colors 
in Figure (3) First, we take a look at the old population. In 
Figure (214., we plot the logarithmic number of stars per ini¬ 
tial mass {\og{dN/dm)) versus initial mass (mini). All stars 
with mini -J 1.0 Mq have evolved to the stellar remnant 
phase, thus they do not contribute to the optical-NIR inte¬ 
grated light any more. The shallower IMF model has smaller 
dN/dm than the standard model in the range of 0.15 to 
1.0 Mq, while the model with 0.3 Mq LMCO has more 
stars in the range of 0.3 to~ 1.0 Mq, due to a normaliza¬ 
tion correction (all stellar population are normalized to the 
same initial mass, so, for example, raising the LMCO means 
that more mass and more light is distributed to the stars 
that remain). The 80% AGE model is not discernible here, 
because the AGE stars have very small initial mass range 
(~ 1.0 Mq). To inspect the strong variance of the post main 
sequence phases, we use the evolutionary points as x vari¬ 
ables instead of mini ■ The evolutionary points are defined as 
points with the same interval arc length along the isochrones 
in the HR diagram. Figure|3)3 clearly shows different number 
of stars for the standard model and the 80% AGE strength 
model. In Figure jSjL, we find that the post main sequence 
stars, though small in number, have major contribution in 
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Figure 1. Color-color and index-index plots for an old, metal-rich population. The IMF slope, LMCO, AGB effects are labelled as 
vectors 1, 2, 3, respectively. At the bottom right of each plot involving photometric colors, an extinction vector of Av^ = 0.02 mag is 
sketched. 
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Figure 2. Color-color and index-index plots for a young, metal-rich population. The IMF slope, LMCO, AGB effects are labelled as 
vectors 1, 2, 3, respectively. At the bottom right of each color-color plot, an extinction vector of Ky = 0.02 mag is sketched. 
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Table 1. Colors and Indices of Each Phase for the Old Population 



PHASE a 

OLD POPULATION 
PHASE b 

PHASE c 

PHASE d 

(B-V) 

1.46939 

0.69314 

0.93189 

1.48725 

(V-K) 

4.04228 

1.75775 

2.87041 

4.84277 

[MgFel 

3.77325 

2.34650 

2.83331 

4.98660 

bTiO 

0.17246 

-0.00274 

0.00573 

0.18546 

Ti02 

0.44303 

0.01672 

0.07097 

0.33609 

WFB A 

2.07944 

0.06787 

0.35408 

1.16696 


the V band integrated light. Then we plot the (B — V), and 
{V — K) colors in Figure[3p and[3}i). Obviously, the very low 
mass stars and the giants have comparable {B — V) colors, 
but the giants, on average, show redder {V — K) colors. 

To better quantify the contribution of different phases, 
we first define four phases: phase a, 0.15 < mini < 0.30 Mq; 
phase b, mini = 0.30 Mq to the main sequence turn-off 
(MSTO); phase c, from MSTO to the beginning of the AGB 
phase, or early AGB (EAGB); phase d, the AGB phase. 
Next, we run the models with only one phase on by nulli¬ 
fying the contribution of all other phases. For example, to 
study the integrated light from phase a stars, we set the 
number of stars in phase b, c, and d to zero. The colors and 
indices from each phase are tabulated in Table [T] A vector 
that connects the standard model and one of these indices 
(or colors) indicates the drift direction if more stars in the 
corresponding phase are added. In Figure |31 these vectors 
are shrunk to one tenth of the original magnitude to accom¬ 
modate the three drift vectors. In the {V — K) vs. {B — V) 
plot, the phase a vector is in the opposite direction of the 
LMCO drift, and similarly the phase d vector has an angle 
of almost 180 degrees with the AGB drift. These indicate 
the LMCO effect is mainly caused by decreasing the contri¬ 
bution of phase a stars, and the AGB effect is mainly caused 
by decreasing the contribution of phase d stars. These two 
conclusions might be obvious, but no vector seems directly 
related to the IMF slope drift. The IMF slope effect changes 
the number of stars of all masses, which means the high-to- 
low-mass ratio of each phase is not a constant. Therefore, we 
need to take the flux contribution of each phase into consid¬ 
eration. In Figure[5)il, we plot the integrated optical spectra 
of four phases. The spectra of phase b and c are both lu¬ 
minous, although phase c stars become brighter than phase 
b stars as wavelength exceeds log(A) « 3.65 (« 4467 A). In 
Figure [SJ3, [Sp, and[Sp, we search for the spectral vari¬ 
ation signals. The spectra of each phase are first normalized 
at log(A) = 3.815. Then the spectrum of the shallower IMF 
model is divided by the standard model. The spectral ra¬ 
tios show shallower IMF model tends to blue the spectrum 
of phase b stars, but redden the spectrum of phase c stars. 
Therefore, the bluing and reddening effects compete again 
each other to determine the integrated color drifts, depend¬ 
ing on the flux ratios of phase b stars to phase c stars (b/c) 
at different wavelengths. Let us take the (B — V) vs. (V — K) 
plot as an example. The b/c flux ratio is close to unity at 
the B and V bands, but much smaller than 1 at the K band. 
Thus the {B — V) color is governed almost equally by both 
the phase b stars and phase c stars, while the {V — K) color 
is mainly controlled by the phase c stars that redden the 


spectrum. The [B — V) color drifts blue due to the stronger 
spectral variation of phase b stars (Figure [5)3) compared to 
phase c stars (Figure (Sp). 

In the index-index plots of Figure|4)3, |4p, and|4p, sim¬ 
ilar conclusions as Figure (dp can be drawn for the LMCO 
and AGB drifts, based on the a and d vector directions. How¬ 
ever, the index variations for the IMF slope effect cannot be 
easily determined from Figure (5) To better quantify the re¬ 
sponse of different phases, we show the index-index plots of 
all four phases in Figure jb] We notice the drift magnitudes 
of the phase b stars are much greater than the phase c stars. 
Since the smallest b/c flux ratio is about 0.5 (at WFB, Fig- 
ure[5p), phase b stars dominate the index variation, pushing 
all the indices to lower values. 

Next, we switch to the young population. A major dif¬ 
ference compared to the old population is that stars with 
mass between 1.0 to ~1.66 Mq still exist (FigureO?). Two 
curves of different slopes intersect at ~1.64 Mq, thus the 
shallower IMF model has more AGB stars than the stan¬ 
dard model! Furthermore, FigurejSp shows longer post main 
sequence phases than the old population. In Figure (3) and 
El, we clearly see the MSTO hook features near the curve 
minimum. The hook feature happens when the star tran¬ 
sitions from core burning to shell burning, with a brief 
gravitational-energy-dominated moment. The {V — K) color 
of the RGB tip is smaller than the old population, due to 
higher stellar mass and thus higher surface temperature of 
the RGB tip stars at young age. In Figure [4)3, |4)^, ET, and 
EH, the directions of vector a and d again testify the previous 
statement that the LMGO effect is mainly caused by phase a 
stars, while the AGB effect is mainly caused by phase d stars. 
To constrain the major contributor of the IMF slope effect, 
we plot the spectra and spectral ratios of all four phases in 
the right column of Figure (5) The spectrum of phase a is 
dwarfed by the flux increase of other three phases. Phase c 
is the most significant of all, starting to overcome phase b at 
log (A) « 3.53. As a result, phase c stars contribute most of 
the light in our optical, NIR colors and indices. Obviously, 
the spectral bluing and reddening competition of phase b 
and phase c stars (Figure E G and H) is won by the latter 
one, leading to redder {B — V) and {V — K) colors for shal¬ 
lower IMF slope. Gomparison of the {B — V) color drifts of 
the old and young populations inspires us to search for the 
critical age where {B — V) color changes from drifting blue 
to drifting red. We find that the (B — V) colors of age less 
than 6 Gyr drift blue, but the {B — V) colors of age greater 
than 8 Gyr drift red. Therefore, between 6 Gyr and 8 Gyr, 
the main contributor of the IMF slope effect for {B — V) 
colors switches from phase b stars to phase c stars. 
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Figure 5. The upper four panels show the flux ratio between the shallower IMF models and the standard models, for four different 
phases. The very bottom panels plot the spectra of the standard models, subdivided into four phases. 


Before ending this section, we calculate the indices of 
each phase (Figure [7]) to explain the IMF slope effects. The 
drift magnitude ratios of phase c over phase b stars has in¬ 
creased to about four times as the old population. Given 
that phase c stars dominate the flux contribution in most 
wavelength, bTiO, Ti02, WFB A, and [MgFe] indices in¬ 
crease. 


3 ADF-IMF COUPLING 

To explore the ADF-IMF coupling, we built a toy model 
in which a hypothetical linear IMF—[M/H] relation is as¬ 
sumed. This linear relation is set by two points: a = 2.35 
at log [M/H] = 0, and a = 1.35 at log [M/H] = —1. Simi¬ 
lar to Paper I, we employ CSPs with single-burst ages but 
normal-width composite abundance distribution functions 
(ADFs). Our normal-width ADF matches well the average 
solar neighborhood ADFs and Milky Way bulge ADFs in 
Paper I. To complete the models, the ADF-weighted Sin¬ 
gle Stellar Populations (SSPs) are combined to construct 
CSPs, in which the IMF slope of each mass bin matches the 
assumed IMF—[M/H] relation. To isolate the consequences 
of assuming a IMF—[M/H] relation, we also construct an¬ 
other CSP models with the normal-width ADF, but instead 
of assuming an IMF—[M/H] relation, we set the IMF to 
be always Salpeter (a = 2.35). We name the latter models 
Constant-IMF CSPs (CCSPs), and call the former models 
Variable-IMF CSPs (VCSPs). To verify the robustness of 
our models, we assemble observed elliptical galaxy spectra 


from Graves et al.l (l2007l ~) , iTraeer et al.l ll2008ll , and IServenI 
||201C ). Readers are referred to Paper I for sample descrip¬ 
tions. Note that all the model and observed indices are cor¬ 
rected to 300 km s”*^ resolution. 

First, we take a look at the optical-NIR color-color plots 
(Figure 1114.). The CCSPs and VCSPs are similar, but show 
differences, especially for the metal-poor populations. Since 
the photometric observations are not included in the above 
samples , we re triev e the (B — V) and ( V — K) colors from 
iPeletierl lll989ll and iPersson et all (Il979ll . 

(i) IPeletieil (|l989l l: In order to avoid focusing on only 
the cores of elliptical galaxies, we choose the colors at maxi¬ 
mum isophotal radius. These colors are corrected for Galac¬ 
tic ex tinction before p lot ting. 

(ii) IPersson et al.l (Il979l l: We choose the well- 
tabulated (B — V) and {V — K) colors of field galaxies. These 
colors have been corrected for Galactic extinction in that pa¬ 
per. 

In Fig.|8}4 the two photometric samples are divided into four 
bins in velocity dispersion (<150, 150—200, 200—250, >250 
km s“^) and we show the median colors of each bin as solid 
stars. These medians confirm that massive elliptical galaxies 
tend have redder {V — K) and [B — V) colors. Note that the 
iy — K) colors are only partially covered by the model grids. 
This model-observation mismatch might be partly because 
CSPs with peak [M/H]> 0.4 is re quired, but we note that 
models with different isochrones llWorthevI Il994l i have no 
trouble reaching red enough {V — K) colors. This is due 
to a slightly stronger coupling between metallicity and the 
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phase a 


age=10Gyr, log[Fe/H]=0.37 
phase b 


phase c 


phase d 
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[MgFe] 
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Figure 6. IMF slope drifts of each phase at old age. 
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Figure 7. IMF slope drifts of each phase at young age. 


temperature of the first-ascent giant branches in the older 
models. The (V — K) miss is almost certainly a model defect. 
We note that that flaw has little bearing on the results of 
this paper which are differential in character. 

Figure |8)3 is the H/f—[MgFe] plot, which is an age- 
metallicity diagnostic diagram. We apply nebular emission 
corrections for H/3 to the Graves sample and Serven sam¬ 
ple following the recipe of IServen fc WorthevI (l2010l l. Five 
emission-corrected galaxies in the Serven sample are labelled 


as open triangles, since larger uncertainties might be ex¬ 
pected for these galaxies in plots involving Balmer features. 
Generally speaking, the observational data points from the 
three samples locate at the old, metal-rich region, as what 
we expected. 

Next, we pick several IMF-sensitive indices to study 
the model variations. Though the Serven sample and the 
Trager sample show moderate scatter compared with the 
model grids, it is encouraging to find the observed data from 


© RAS. MNRAS OOO.ITIflil 






























8 Baitian Tang and Guy Worthey 


age=10 Gyr, log[Fe/H]=0.37 age=2 Gyr, log[Fe/H]=0.37 






Evolutionary points Evolutionary points 



Evolutionary points 



Evolutionary points 




Evolutionary points Evolutionary points 




Evolutionary points Evolutionary points 


age=10 Gyr, log[Fe/H]=0.37 


age=2 Gyr, log[Fe/H]=0.37 





[MgFe] 



[MgFe] 



[MgFe] 



[MgFe] 


Figure 4. To illustrate the connections of effects and evolutionary 
phases, we connect the standard models with the contributions of 
each phase. The phase vectors (a,b,c,d) are shrunk to one tenth of 
the original magnitude to accommodate the three effect vectors. 
The vector of each phase is labelled by the corresponding phase 
name. 


Figure 3. The number of stars, luminosities, and colors of the 
old, and young populations are plotted against isochrone evolu¬ 
tionary points or mini , where the use of evolutionary points better 
higlights post main-sequence phases. Four phases are defined as: 
phase a, 0.15 < mini < 0.30 Mq; phase b, mini = 0.30 Mq to 
the main sequence turn-off (MSTO); phase c, from MSTO to the 
beginning of the AGB phase, or early AGB (EAGB); phase d, the 
AGB phase. 


[Graves et alj (|2007l) are well-covered by our models, and the 
variation trend looks reasonable. 

In all the panels of Figure [S] we find that the VCSPs 
appear more metal-rich than the CCSPs. According to our 
IMF—[M/H] relation, the IMFs of the metal-poor popula¬ 
tions in the VCSPs are shallower. We show that the shal¬ 
lower IMF model appear dimmer than the standard IMF 
model in Figure [S] Therefore, the metal-poor populations 
appear dimmer in the VCSPs. Note that metal-poor popu¬ 


lations outshine their metal-rich counterpart^, and thus the 
smaller integrated-light contribution from the metal-poor 
populations leads to more metal-rich appearing VCSPs. 


4 DISCUSSION 
4.1 Swapping Models 

To test the robustness of the model drifts, we also per¬ 
form similar experime nts with the Flexible Stellar Popula¬ 
tion Synthesis fFSPS. IConrov &i Gunnllioioh . We arranged 
to utilize the FSPS with the same age, IMF slope, and 
LMCO as our models. But the maximum metallicity of the 
Padova-I-MILES option is log[Fe/H]= 0.20, which is differ¬ 
ent than the metallicity we chose above (0.37). All the IMF- 
related spectral changes might be age and metallicity sensi¬ 
tive. Since the wavelength coverage of MILES is 3600—7400 
A, the Nai and WFB A are not available for analysis. To 

® If both populations have similar number of stars 
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[MgFe] 
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[MgFe] 



[MgFe] 



1 - IMF slope 

2-LMCO V 

0.072 


3-AGB 

O 0.070 

0.068 

3 


0.066 

2 


0.064 



3.24 
3.22 
3.20 
3.18 
3.16 
0.960 0.970 


0.980 

B-V 


0.990 1.000 


5.14 5.16 5.18 
[MgFe] 



0.082 



0.080 



CM 0.078 
O 

2^ 

3 

1- 0.076 



0.074 



[MgFe] 



5.14 5.16 5.18 
[MgFe] 


Figure 8. The VCSPs (dashed lines) and CCSPs (solid lines) 
are described i n [3] In the top left panel, me asuremen ts from 
IPersson et alJ lll979l 'l (open pentagons), and IPeletiei] lligHfli l 
(crosses) are shown. The medians of different velocity dispersions 
(<150, 150—200, 200—250, >250 km s“^) are indicated as solid 
stars. In the optical index plots, Serven (triangles). Graves (open 
diamonds), and Trager (squares) observations are shown. Five 
Balmer emission galaxies in the Serven sample (open triangles) 
have had additional corrections applied. 


evaluate the AGB effect in FSPS, we choose the o riginal 
Padova TP-AGB treatment llMarigo fc Girar^l2007l ) as the 
standar d model, and compare it with the TP-AGB treat¬ 
ment of IConrov fc Gunnl J201fll ). Note that the latter TP- 
AGB treatment effectively reduces the NIR related colors 
at age around 1 Gyr, and well fits the colors of star clus¬ 
ters in the Magellanic Cloud (See their Figure 3). But both 
TP-AGB treatments have similar descriptions at old ages, 
and thus we expect small AGB drift for the old population 
in FSPS. This is different than the way we change the AGB 
strength in 112.21 the AGB strength is always reduced to 80% 
of the full strength, regardless of the age. 

We show the color-color plots and optical index-index 
plots in Figure [9] and [To] To compare with the results from 
our models, all the dynamic ranges are set to be the same 
as Figure [T] and [2| It is encouraging to find a lot of similar¬ 
ities in the directions and magnitudes of the drifts between 
two sets of models, but discrepancies still exist. For the old 
population, all the AGB drifts are very small, as expected. 
The IMF slope drifts seem smaller in the (F—K) colors, but 
show much bigger magnitudes in the negative [MgFe] direc¬ 
tion. Note that the LMCO drift is no longer greater than 
the IMF slope drift for the Ti02 index: the IMF slope effect 
is the leading effect, now. Turning to the young population 
model, the color drifts of LMCO and IMF slope seem ro¬ 
bust between two sets of models, but the AGB drifts are 
greater in the FSPS models, due to the strong reduction of 
TP-AGB stars around 1 Gyr. Surprisingly, the index drifts 


Figure 9. Color-color and index-index plots at old age using 
FSPS models. 
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Figure 10. Color-color and index-index plots at young age using 
FSPS models. 

are much smaller than our models. This is possibly because 
of the different index derivation routines of these two model 
sets, but we are not sure of that. 

4.2 Feasibility of Breaking the Degeneracy 

Exploring whether the apparent steep IMF in massive ellip¬ 
tical galaxies might be partially due instead to the decreas¬ 
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ing number of AGB stars requires one to be able to find an 
observational signature distinguish the two effects. We add 
to that list the LMCO, of course, but also note that various 
age and metallicity effects need to also be addressed, his- 
torically derived from the H/3— <Fe>PI or Hyrr— <Fe> plot 
dTrager et al.ll200^ y: iTang et ^l2009l ~). 

For a young elliptical galaxy, the color-color plots of Fig- 
ure[^show substantial drifts (~0.03 mag) for the IMF slope 
and AGB effects. Note that the AGB and IMF slope effects 
vector oppositely in most of the Figure [2] plots . Since our 
results should be interpreted in a partial derivative sense, 
the opposite vectors in fact do nothing to disentangle the 
IMF slope and AGB effects. That is to say, the degeneracy 
is still firm and hard to break in this case. 

But for an old, metal-rich galaxy, the color-color plots 
of Figure [T] show the LMGO, IMF slope, and AGB effects 
are distinguishable if the photometric accuracy is better 
than 0.02 mag. The index-index plots reveal robust IMF 
slope drift in the [MgFe] direction, reaching an amplitude of 
0.02 A. This separates the IMF slope effect from the LMGO 
effect, for the latter one causes much smaller drifts in the 
[MgFe] direction. Note that the [MgFe] index is insensitive 
to [a/Fe], giving the [MgFe]-related plots advantage to dif¬ 
ferent alpha-enhanced environment. 

We find therefore that it is practical to lift the degen¬ 
eracies using optical-red spectra and optical-NIR photome¬ 
try if the measurements are accurate enough. We estimate 
these accuracies should be achieved: (1) dust attenuation 
o'Ay <0.01 mag; (2) photometric accuracy a <0.01 mag; 
(3) index accuracy a <0.02 A. As s tellar photometric ac¬ 
curacy reaches milli-mag these days llClem et al.ll200u L an 
accuracy of 0.01 mag is feasible for nearby galaxies if one 
is precise about the observational complications. To list a 
few of them: dust attenuation, both Galactic and in-situ; 
sky subtraction for extended sources; aperture matching for 
spectroscopic and photometric extractions. 

We also notice that the Nai index seems to be insensi¬ 
tive to the AGB effect for both the old and the young pop¬ 
ulations, in spite of Nal’s red central wavelength. A look 
at the index fitting functions at the stellar level reveals this 
insensitivity may relate to the similar Nal indices of giants 
and dwarfs. The Nal index might be a way out of the ADF- 
AGB-IMF masquerading. 


4.3 Symposium of Multiple Effects 

Age and metallicity effects are generally strong, but also 
degen erate in most of the colors and indices (IWorthevI 
I 1994 I I. To obtain a bigger picture besides the three ef¬ 
fects that we discuss in this paper, we summarize the ef¬ 
fects concerning the IMF slope, LMGO, AGB contribution, 
age, metallicity, ADF, and VGSP in Figure 1111 They are 
labelled as arrow 1 to 7, respectively. The parametriza- 
tion of the first three effects is described in 1)2.11 The age 
effect is estimated by calculating SP models of age=12 
Gyr, log[Fe/H]=0.37, and age=4 Gyr, log[Fe/H[=0.37 (A 
age=-|-2 Gyr), while the metallicity effect is estimated by 
calculating SP models of age=10 Gyr, log[Fe/H]=0.185, and 


<Fe>=(Fe5270 -I- Fe5335)/2 llGonzalejll^ . 


age=10 Gyr, log[Fe/H]=0.37 age=2 Gyr, log[Fe/H]=0.37 








Figure 11. To visualize effects that may impact the SSP param¬ 
eter determination, IMF slope, LMCO, AGB contribution, age, 
metallicity, ADF, and VCSP effects are labelled as vector 1 to 7, 
respectively. 


age=2 Gyr, log[Fe/H]=0.185 (Alog[Fe/H]= —0.185). To in¬ 
vestigate the effects concerning GSPs, we move the models 
with normal-width ADF, peak [M/H]= 0.4, and constant 
IMF slope to the center. The ADF effect is shown as a dis¬ 
placement from the normal-width ADF GSPs to the narrow- 
width ADF GSPs, and the ADF-IMF coupling is represented 
by a vector from the GCSPs to the VGSPs. 

Age and metallicity effects are parallel and degenerate 
as expected. They are joined by two new effects concerning 
the GSPs: the ADF and VGSP effects. In fact, the ADF ef¬ 
fect is the same as the “red lean” effect (Paper I) where nar¬ 
rower ADF appears more metal-rich than a wider one. The 
reasons for the VCSP effect are explained in 331 Careful read¬ 
ers may find that the ADF and VCSP effects both originate 
from the suppression of metal-poor populations, which leads 
to a metal-richer appearing integrated-light model. The sim¬ 
ilar origins and vector directions drive us to segregate the 
age, metallicity, ADF and VCSP effects into one group, and 
label them as group I effects. Group II effects are IMF slope, 
LMCO, and AGB contribution. Group II effects vector al¬ 
most orthogonally to group I effects. Though group I effects 
are more prominent in magnitude, the orthogonality sug¬ 
gests that observations to isolate group II effects are still 
feasible. 
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4.4 Recovering [X/Fe], IMF slope, LMCO, and 
AGB Percentage 

IWorthev et alJ (l2014lj) showed our efforts to recover elemen¬ 
tal abundances from observed spectra. In this work, we up¬ 
grade our inversion program by replacing the SSP models 
with normal-width ADF CSP models and adding the IMF 
slope, LMCO, and AGB percentage (AGB%) parameter de¬ 
termination. 

To examine the nncertainties of our inversion program, 
500 mock galaxy spectra are constructed by Monte Carlo 
simulation. First, a model is produced with known age, 
abundance, IMF slope, LMCO, and AGB fraction. In this 
model the indices and colors are altered by random numbers 
which have a normal distribution. The standard deviation 
of the normal distribution function is set to O.Lustd, where 
(Tstd is a vector composed for reasonable guesses for the ob¬ 
servational uncertainty in each color or index (Tabled]). We 
supply the inversion program with each of these mock galaxy 
spectra in turn, and compare the recovered values with in¬ 
put values. Figure [T^l show the recovered and input values 
of models at age = 12 Gyr, peak [M/H] = —0.1. 

In FigurefT2K. the age and peak [M/H] are quantized in 
the current inversion algorithm, since these two basic param¬ 
eters are determined by closest-match rather than a smooth 
interpolation within the model grid. Age-metallicity degen¬ 
eracy is the reason for the anti-correlation found in the first 
panel. In Figure fT2l 3. we see tight correlation between [C/R] 
and [N/R], since these two elements are largely determined 
by common indices, the GN bands around 4100 A. In Fig¬ 
ure [HP, HUT, and HUH, the correlations among IMF slope, 
LMCO, and AGB% are generally loose. This implies the de¬ 
generacies among these three parameters are, for the most 
part, broken by our selection of indices. Though modest scat¬ 
ter exist in the recovered values, all of the mean recovered 
values agree with the input values inside the error range 
without alarming systematics. Some of the mean recovered 
values are especially close to the input values, e.g., IMF 
slope, [O/R] and [Na/R]. This success encourages us to ap¬ 
ply our inversion program to observed galaxy spectra and 
colors in the future. 


5 SUMMARY 

A modeling study of various new parametric effects on in¬ 
tegrated light, begun in Paper I by exploring the effects 
of the Abundance Distribution Function (ADF), contin¬ 
ues on in this paper to explore two underexplored effects 
that might mimic IMF slope variation: the LMCO and 
the AGB strength. If we hypothesize that the apparently- 
steep IMF in massive elliptical galaxies might be partially 
due to the decreasing number of AGB stars, the ques¬ 
tion is can we ever know it, given the stellar population 
degeneracies. The magnitudes of the IMF slope, LMCO, 
and AGB effects are smaller than age and metallicity ef¬ 
fects, but these two groups of effects vector almost orthog¬ 
onally in diagnostic diagrams, which bodes well for mea¬ 
suring the combination of slope/LMCO/AGB effects. Inter¬ 
nal to that trio, we explore the degeneracies and find that 
it is very difficult to distinguish steepening IMF from de¬ 
creasing AGB strength for young, metal-rich populations. 


However, the slope/LMCO/AGB degeneracies can be lifted 
for old (age«10 Gyr), metal-rich populations using a com¬ 
bination of optical-to-near infrared photometry and spec¬ 
troscopy. The disentanglement happens, however, at an ob- 
servationally challenging level («0.02 mag). 

We fragmented our models into different evolutionary 
phases to isolate and discuss the leading factors for the IMF- 
related effects. We also investigated a series of models with 
an ADF-IMF coupling in which metal-rich populations fa¬ 
vor low mass star formation. Models with ADF-IMF cou¬ 
pling appear more metal-rich than the noncoupled models, 
strongly resembling a modest narrowing of the ADF width, 
which has the same “red lean” effect. Attempting to recover 
the magnitude of an ADF-IMF coupling is very challenging 
from integrated light alone. _ 

W e upgraded our inversion program llWorthev et al.l 
l2014bfl by replacing the SSP models with normal-width ADF 
GSP models and adding the IMF slope, LMCO, and AGB 
percentage parameter determination. We estimated uncer¬ 
tainties on parameter estimation from Monte Carlo simu¬ 
lations, and find no significant systematic drifts, though of 
course we see many partial parameter degeneracies. The de¬ 
generacies among the IMF-related parameters can be broken 
by our selection of indices for old stellar populations. This 
success encourages us to seek appropriate observational ma¬ 
terial to apply our techniques to observed galaxy spectra 
and colors in the future. 
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Table 2. Baseline Assumed Observational Uncertainty in Each Color or Index 



Index 

(7 

Unit 


Index 

<7 

Unit 


Index 

cr 

Unit 


Index 

(7 

Unit 

1 

CNl 

0.005 

mag 

24 


0.121 

A 

47 

Cr3594 

0.131 

A 

70 

bTiOsTXC,^ 

0.005 

mag 

2 

CN2 

0.006 

mag 

25 

H7f 

0.114 

A 

48 

Cr4264 

0.157 

A 

71 

aTiOsTKC 

0.003 

mag 

3 
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0.081 

A 

26 

Ho 

0.114 

A 

49 
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A 

72 
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0.004 

mag 
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A 
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0.339 

A 

50 
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A 

73 
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mag 
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A 
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A 

51 
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A 

74 
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A 
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A 

29 
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0.221 

A 

52 
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A 

75 

Ti02sDSst, 

0.003 

mag 
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Fe4531 
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A 

30 

CN04175 

0.285 

A 

53 

Mn4757 

0.135 

A 

76 

H/3emiss^ 

0.090 

A 
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0.319 

A 

31 
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A 

54 

Fe4058 

0.173 

A 

77 
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A 
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A 
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A 
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A 
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mag 
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mag 
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mag 
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mag 
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mag 
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A 
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A 
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0.113 

A 

91 

(U-K) 

0.020 

mag 

23 

H7a 
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A 
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A 

69 
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Figure 12. The recovered values (scattered dots) are compared with the input values (solid lines). The mean recovered values and their 
corresponding standard deviations are labelled with error bars. 
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